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SUMMARY 


A supersonic tunnel designed to operate at Mach number 1.3 was used 
for a preliminary experimental flutter investigation of widely different 
unswept cantilever wings. Data for 12 wings with mass-density parameters 
1 /k ranging from 52 to 268 , center -of -gravity positions ranging from 4-6 
to 63 percent chord from the leading edge, and elastic -axis positions 
ranging from Jh to 52 percent chord from the leading edge are considered. 

A comparison is made of the test results with calculations of 
bending-torsion flutter obtained by the theory of flutter in supersonic 
two-dimensional flow and it is concluded that the test data are in reason- 
able agreement with the calculated results. In general, the theoretical 
values are conservative. As shown by the theory, the flutter results are 
quite sensitive to the location of the center of gravity. Thick and thin, 
blunt and sharp airfoil-section shapes were used, but no very pronounced 
effect of the section shape on flutter characteristics was found. The 
experiments include a study of the effect of the addition of tip moments 
of inertia. With the center of gravity of the tip weights coincident 
with the center of gravity of the wing section, no detrimental effect on 
the flutter speed was found. 


INTRODUCTION 


The background and theory for the flutter of an airfoil in a two- 
dimensional flow at supersonic speeds is given in reference 1. The pres- 
ent investigation is a preliminary survey to determine the possibility 
of using the theory of reference 1 for flutter at supersonic speeds to 
predict the coupled bending-torsion flutter of widely different unswept 


* Supersedes the recently declassified NACA RM iBJll, "initial ^ 
Experiments on Flutter of Unswept Cantilever Wings at Mach Number 1.3 by 
W. J. Tuovila, John E. Baker, and Arthur A. Regier, 19^9* 
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cantilever wings at a low supersonic Mach number. This preliminary inves- 
tigation is not intended as a critical test of the theory since the analy- 
sis does not consider the effect of mode shape, aspect ratio, section 
shape, tip Mach cone, or viscous effects. 

A single-degree-of -freedom torsional instability which may occur in 
the Mach number range 1.0 to I .58 is discussed in reference 1. In order 
also to investigate the possible occurrence of such single-degree flutter 
on cantilever wings, the test apparatus was designed to operate at a Mach 
number of 1.3. 

A brief discussion is given of the effects of concentrated masses 
placed at the wing tip, the center of gravity of the masses coinciding 
with the center of gravity of the wing, and the effects of sharp and 
blunt leading edges on the wings. 


SYMBOLS 

b semichord, ft 

c w chord, ft 

GJ torsional stiffness, in2-lb 

p mass density of air in test section, slugs/cu ft 

m mass of wing per unit span, slugs/ft 

M Mach number 

1/k mass-density parameter, m/npb^ 

I a mass moment of inertia of wing about elastic axis per unit span 

r a 2 radius -of -gyration parameter, Ia/mb 2 

V flutter velocity, ft/sec 

ff flutter frequency, cps 

fh first-bending frequency, cps 

f a first-torsion frequency, cps 


(Df* = 2rtff 
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= 2itf h 

= 2jtfa 

gk first-bending damping coefficient 

first-torsion damping coefficient 


APPARATUS AND TEST METHODS 


The tests were made at a Mach number of 1.3 in an M intermittent” 
two-dimensional supersonic tunnel having a 9*24-inch by 18. 23-inch test 
section (figs. 1 and 2). The apparatus operates from atmospheric pressure 
to a vacuum. A quick-operating valve allows a steady-flow condition to 
be reached in approximately 0.15 second and this condition of steady flow 
can be maintained for 5 to 7 seconds. To prevent condensation in the test 
section, a room was constructed at the tunnel entrance in which the air 
could be heated. Variation of the air temperature from I7O 0 F to 200° F 
results in a velocity range at the test section from 950 miles per hour 
to 990 miles per hour at a Mach number of 1.3- The test-section Mach num- 
ber determined by optical means varied from 1.29 to 1.31. The test- 
section Mach number determined by a pressure survey showed a variation 
from 1.31 to 1.33 (fig* 3)* 

The models were mounted cantilever fashion from the side of the 
tunnel (fig. 4). In order to avoid possible model failure during the 
tunnel transient conditions, and since supersonic flutter characteristics 
were being determined, it was desirable to withhold the model from the 
flow during the periods of acceleration and deceleration. A pneumatic - 
cylinder arrangement was installed for this purpose (fig. 2). With this 
device, the model could be held outside the tunnel wall until stable flow 
conditions were reached; then, by means of electrically operated valves 
controlled by an observer, the model could be projected into the tunnel 
slowly and withdrawn quickly if necessary. 

The flutter models were of rectangular plan form and were constructed 
of laminated Sitka spruce or aluminum alloy. The wing dimensions and 
parameters are given in table I. The wing chords ranged from 3-03 to 
4.22 inches and the lengths (semispan) from 6 to 9*12 inches. Both thick 
and thin sections with blunt and sharp leading edges were used. The air- 
foil sections used were 3-, 5-> and 8-percent-thick circular arcs, 

3 -percent-thick double wedge, NACA 16-010, and NACA 65-OO7. The mass- 
density parameter 1 /k ranged from 52 to 268, the chordwise positions of 
the center of gravity ranged from 46 to 63 percent chord, and the posi- 
tions of the elastic axis ranged from 34 to 52 percent chord. 
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Before each model was installed in the tunnel, its weight, moment of 
inertia, and section center-of -gravity position were determined. After 
installation in the tunnel, the elastic axis was located by observing the 
axis of zero twist optically. The first-bending frequency and the damping 
were obtained from a free-vibration record of the model. Since the wings 
were structurally uniform, the stiffness parameter GJ could be computed 
from a torsional-vibration record obtained with a mass of large known 
moment of inertia attached to the wing tip. The uncoupled first -torsion 
frequency could then be calculated. The torsional damping was determined 
from the torsional-vibration-decay curve. 

The models were designed not to flutter on the first run. The models 
were progressively modified by adding lead to shift their centers of 
gravity and by slotting to shift their elastic sixes and change their 
bending and torsion frequencies, 'until flutter occurred. If a model was 
saved, its parameters were changed to yield another flutter point. The 
values of the radius -of -gyration parameter ra 2 and chordwise center of 
gravity were determined from the unmodified wing and the added masses. 

During each test run, the following data were recorded simultaneously 
by means of a recording oscillograph: flutter frequency, position of the 

model in the tunnel, and static pressure. 

A sample record of the flutter of model B-5 is given in figure 5- 


METHOD OF ANALYSIS 


The flutter data for the particular wings tested are put in coeffi- 
cient form and compared with the analytic solution of the supersonic 
bending-torsion flutter problem in a two-dimensional flow given in ref- 
erence 1. The data of this paper were obtained at a Mach number of 1.3 
and, since aerodynamic coefficients at this Mach number are not included 
in reference 1, a linear interpolation was made between the computed 
values of the flutter-speed coefficient at Mach numbers of 1.25 and 1-43, 
for which the aerodynamic coefficients are available. First-bending and 
uncoupled first-torsion frequencies and damping coefficients gh and 
were used in the theoretical analyses. 

An examination of the possible errors introduced into the results by 
the method of interpolation is desirable. It is known that the error may 
be very large; for example, in the case of torsional instability in one 
degree of freedom for the elastic -axis position at 50 percent chord, the 
interpolation was made directly between the aerodynamic coefficients at 
Mach numbers 1.25 and 1.43. This interpolation was necessary since the 
wing was stable at a Mach number of 1.43 and the theory yields no solution 
for the flutter-speed coefficient under these conditions. 
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RESULTS AND DISCUSSION 


The significant flutter parameters and the results of the calcula- 
tions are given in table I. A large number of tests were made on wings 
which did not flutter, but table I gives only the results for the wings 
for which flutter was observed. Altering a model to obtain flutter con- 
sisted in moving the center of gravity back in increments of about 2 per- 
cent of the chord. Consequently, the chordwise position of the center 
of gravity that would produce flutter is known to within about 2 percent 
chord. The results are sensitive to center -of -gravity position and this 
sensitivity may account for some of the scatter of the data. Contributing 
also to the scatter of the data are the inaccuracies in obtaining the wing 
parameters, the effect of the degree of penetration of the model into the 
tunnel, and errors in the determination of the flutter-speed coefficients 
arising from the method of interpolation. 

Some of the results listed in table I are presented in figures 6 
and 7. In figure 6 the theoretical and experimental flutter-speed coef- 
ficients are compared. The fact that the experimental data fall above 
the 4 5° line, for the most part, indicates that the theory of reference 1 
is generally conservative as far as application to cantilever wings is 
concerned. From consideration of the fact that a slight inaccuracy in 
the location of the center of gravity has a large effect and also that 
effects of section shape, aspect ratio, and Mach cone are not accounted 
for, the agreement is not unsatisfactory. The theoretical and experi- 
mental flutter frequencies are compared in figure 7} the experimental 
frequencies ranged from about 0.6l to 1.09 times the theoretical values. 

In all cases the mode at flutter appeared to consist mainly of a coupling 
of first-bending and first-torsion modes. 

Since the present investigation is of a preliminary nature and covers 
a wide range of parameters, no attempt was made to isolate the effects of 
separate parameters such as the mass -density parameter, frequency ratio, 
elastic axis, and center of gravity, which are treated by the two- 
dimensional theory, or parameters such as aspect ratio not treated by 
the theory. 

An attempt was made to investigate some of the possible effects of 
airfoil-section shape on flutter. The intermingling of the data for the 
various models (figs. 6 and 7) suggests that the section shape has no 
very pronounced effect on flutter at a Mach number of 1.3. However, more 
difficulty due to divergence was experienced with thick models having 
blunt leading edges. This result is in accord with the higher-order 
method of calculation (order higher than in the linear method) for pres- 
sure distribution at supersonic speeds in steady flow, which shows that 
the center of pressure may move ahead of the 50-percent-chord position 
for thick blunt-nosed airfoils, particularly at Mach numbers near unity. 
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It was observed in the tests that the thick airfoils tended to diverge 
even though the elastic axis was ahead of the ^“Percent-chord position. 

Since practical winged vehicles pass through the subsonic speed range 
in order to reach supersonic speeds, some discussion of and comparison 
with subsonic criteria are desirable. For this purpose, incompressible 
flutter-speed coefficients were computed by the method of reference 2 
for which first bending and uncoupled first torsion frequencies and 
damping coefficients gh and go, were utilized. That flutter-speed 
coefficients based on subsonic theory are conservative with respect to 
supersonic speeds has been suggested in reference 3 and, in fact, is 
indicated by some of the calculations in reference 1. This conclusion 
is also indicated by the present tests, as illustrated in figure 8, in 
which the experimental flutter-speed coefficients are plotted against 
the incompressible theoretical flutter-speed coefficients. The statement 
may not be true in general; for example, the condition when the frequency 

ratio ^ = 1 may need further investigation and, in any case, the margin 

of safety is not established. Some of the models were permitted to 
encounter the tunnel transient speeds and, for example, model F-l, which 
had fluttered at Mach number 1.3> was held in the tunnel while the tunnel 
was brought up to speed. The wing fluttered and broke at a Mach number 
of about 0.7, a result which is in general agreement with the subsonic 
calculation. Flutter data obtained with bombs and rocket missiles (ref- 
erences 4 to 6) and other experience indicate that if flutter failures 
occur, they occur, in general, at speeds near sonic. For the practical 
purpose of making preliminary estimates of a wing flutter speed, such 
formulas as, for example, the approximate flutter formula in reference 2 
or similar criteria thus appear useful over a wide range of speeds. 

In reference 3, Smilg suggests a torsional frequency criterion 
ofc(C w > 1,000 feet per second as sufficient to prevent one degree of tor- 
sional flutter at supersonic speeds. The criterion is based on the 
assumption that the first-bending frequency is very high with respect 
to the first-torsion frequency. In order to explore this criterion, a 
cantilever model was equipped with tip weights at both the leading and 
trailing edges so that the torsional frequency was reduced. The results 
of the tests are shown in table II. In all cases a slight shift of the 
center of gravity ahead of the location at flutter stopped the flutter. 

The fact that flutter is extremely sensitive to the center -of -gravity 
position and that the values of the product c%c w are far below 1,000 
indicates that the criterion is overly conservative when applied to 
cantilever wings with normal bending-torsion frequency ratios. The data 
further suggest that for cantilever wings the bending degree of freedom 
may suppress the one-degree-of -freedom torsional flutter and that bending- 
torsion effects occur. 
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The data further Indicate that no harmful effect on the flutter 
speed ensues when the center of gravity of the tip weights and the wing 
coincide. It is observed that the frequency ratio varies from 0.55 
to 0.85 and that the torsional frequency has been reduced to as low as 
one-third of the value for the wing without tip weights. For the largest 
mass moment of inertia on the wing tip (last case in table II), it was 
necessary to move the center of gravity farther toward the trailing edge 
to produce flutter. 

In figure 9> the theoretical curves represent flutter-speed coeffi- 
cients for one -degree -of -freedom flutter calculated according to refer- 
ence 1. The experimental flutter-speed coefficients shown in this figure, 
however, correspond to the coupled bending -tors ion values. 

An effort was made to obtain some systematic aspect-ratio effects 
from the present tests, but the results were rather contradictory. Some 

models which spanned the tunnel (except for — — inch tip clearance) were 
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used so that, presumably, two-dimensional flow over the wing could be 
expected. Flutter of full-span models of NACA 16-010 section could be 
stopped by retracting the tip from the boundary layer; however, when the 
tip was retracted from the boundary layer for the 3 -percent -thick double- 
wedge models, the flutter amplitude increased. The effect of the subsonic 
boundary layer at the tip of the models is not known. In one particular 
case the model fluttered at 99 percent of the theoretical frequency on 
entering the tunnel and the frequency changed to 128 percent at a smaller 
amplitude when the model spanned the tunnel. As the model was retracted, 
the flutter frequency dropped to 99 percent of the theoretical value and 
fluttered to destruction. A more systematic investigation of the aspect 
ratio and tip and shape effects is desirable to clarify various aspects 
of the problem. 


CONCLUSIONS 

The results of a preliminary flutter investigation of widely differ- 
ent unswept cantilever wing models at a Mach number of 1.3 indicated the 
following conclusions: 

1. Agreement between experimental and calculated flutter-speed coef- 
ficients is fairly satisfactory. In general, the theoretical values are 
conservative. 
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2. No very pronounced effect of airfoil-section shape on the flutter 
characteristics was found in these preliminary experiments; however, sig- 
nificant divergence effects were observed on thick blunt-nosed airfoils. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 13, 19^-8. 
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TABLE I 


EXPERIMENTAL AND THEORETICAL RESULTS OF FLUTTER INVESTIGATIONS 


— ~ -____lfodel designation 

Parameters " — 

A— 1 

B— 1 

B— 2 

B— 3 

B— 4 

B-5 

C-l 

C— 2 

D— 1 

E— 1 

F— 1 

G— 1 


NACA 

NACA 

NACA 

NACA 

NACA 

NACA 

8 percent 

8 percent 

Modified 
circular arc 

3 percent 

3 percent 

Modified 
circular arc 
5 percent 

Description of model section . . . 

65-007 

16-010 

16-010 

16-010 

16-010 

16-010 

circular arc 

circular arc 

4.74 percent thick 

circular arc 

double wedge 

thick 

Length, in 

6.06 

9.125 

9.125 

9.125 

9.125 

7.50 

6.00 

6.00 

7.125 

9.125 

9.125 

6.00 

Chord, in 

3.03 

4.03 

4.03 

4.03 

4.03 

4.03 

3.03 

3.03 

4.22 

4.04 

4.01 

4.00 

Center of gravity, percent chord . 

U9.1 

51.6 

54.6 

56.0 

56.7 

57.0 

53-0 

55.4 

46.0 

62.6 

56.5 

53.5 

Elastic axis, percent chord . . . 

Ul .3 

3^-1 

39.6 

39.6 

44.2 

39-5 

48.0 

51.55 

37.0 

38.7 

45.2 

47.5 


0.28 

0.39 

O .38 

0.40 

0.37 

0.37 

0.230 

0.233 

0.275 

0.510 

0.29 

0.27 

1 / K 

64.9 

95.3 

108.1 

113.1 

113.3 

130 

67.1 

74.1 

53.5 

267.5 

150.8 

51.7 

GJ 

578 

3060 

2945 

3150 

4105 

2500 

839 

909.5 

1025 

2710 

2220 

1066 

V 

133.4 

95.0 

100.9 

98 .O 

104.2 

104.0 

157 

157 

no 

25.2 

28.5 

148.5 

V CpS * • * ‘ 

278 

163 

156.3 

154.8 

183 

162 

363 

361 

178 

81.7 

132.8 

245 

f f , cps 

180 

134 

132.5 

134.6 

131.5 

136 

188 

184 

142 

71 

70.5 

176 

p x 10 ^ (test section) 

0.876 

0.945 

0.890 

0.897 

O .908 

0.886 

0.887 

0.884 

O .918 

0.888 

0.888 

0.910 

«h 

0.04 

0.03 

0.025 

0.03 

0.025 

0.035 

0.05 

0.05 

0.03 

0.004 

0.01 

0.04 

«a 

0.04 

0.04 

0.04 

0.035 

0.03 

0.035 

0.04 

0.04 

0.03 

0.005 

0.005 

0.04 

“h/“a ' 

0.48 

0.583 

0.645 

0.633 

0.57 

0.64 

0.432 

0.435 

0.614 

O .308 

0.215 

0.606 

V/bce f (experimental) 

10.15 

9.98 

10.31 

10.20 

10.40 

10.03 

9.715 

9.92 

9.04 

19.13 

19.61 

7.71 

V/baij. (theoretical) 

6.946 

7.770 

7.109 

7-546 

7.938 

7.683 

8.53 

7.18 

5.28 

14.85 

9.38 

4.92 

V/bco a (experimental) 

6.59 

8.21 

8.74 

8.91 

7.44 

8.42 

5.02 

5.055 

7.25 

16.7 

10.35 

5.54 

V/bo) a (theoretical) 

6.554 

7.158 

5.664 

6.029 

6.234 

6.041 

4.913 

5.503 

5.069 

13.08 

8.23 

3.905 

(experimental) 

0.648 

0.822 

0.847 

0.870 

0.719 

0.840 

0.518 

0.51 

0.798 

0.868 

0.531 

O .718 

a (theoretical) 

0.944 

0.921 

0.797 

0.800 

0.785 

0.786 

0.577 

0.766 

0.96 

0.88 

0.868 

0.795 

(V/bu> a ) (M ■ 0) (theoretical) . . . 

3*427 

4.302 

4.245 

4.358 

4.121 

4.612 

3.184 

3*374 

2.811 

8.263 

5.491 

2.585 
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TABLE II 

FLUTTER DATA ON A CANTILEVER MODEL WITH TIP WEIGHTS 

Wing chord = 4 in.; elastic axis at 47 percent chord; r a 2 = 0.23 (for wing with 
no tip weights); weight of wing with no tip weights = 0.0806 lb; length of 
wing (semispan) = 6^ in.; tip-weight center of gravity coincides with the 
wing center of gravity 


Frequency, 

cps 


Center of gravity 
of wing and 
tip weights, 
percent chord 

fps 

Tip weights, 
lb 

Moment of 
inertia of 
tip weights 
about c.g., 
in-lb-sec^ 

Torsion 

Bending 

Flutter 

L. E. 

T. E. 

244 

134 

153 

O.550 

50 

510 

0 

0 

0 

133 

102 

105 

.767 

50 

278 

.00949 

.00949 

.000162 

103 

81 

86 

.786 

49 

216 

. OI766 

.01754 

.000342 

80 

68 

74 

00 

VJ1 

O 

53 

167 

.02747 

.03044 

.000686 
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Figure 1.- Diagram of the supersonic flutter apparatus. 
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Figure 2 . - General view of the Langley supersonic flutter apparatus . 
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Figure 3 • - Mach number 3urvey of the test section. 
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V/bo^ theor 


Figure 6.- Comparison of experimental values of V/bo^ 

values of V/bc%, at M = 1.3. 


to theoretical 
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Figure 7.- Comparison of experimental values of cuf/o^x 

values of otf/oa at M = 1.3. 


with theoretical 
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Figure 8.- Comparison of experimental values of V at M = 1.3 
with theoretical incompressible values of Y /bcu a at M = 0. 
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Figure 9 • - Comparison of experimental values of V/Tx^ for coupled wing 

flutter with theoretical values for one -degree- of -freedom torsional 
flutter. M = 1 . 3 , g a = 0.04. 
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